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ABSTRACT 

We present a new model for low-luminosity X-ray sources in globular clusters, with Lx < 10^"* erg s~'. The 
model we propose is that of a single neutron star accreting from cluster gas that has accumulated as a natural 
product of stellar evolution. An analytic luminosity function is derived under the assumption that the speed 
distribution of neutron stars in the central region of a cluster is described by a Maxwellian, and that the density 
and temperature of the gas are uniform. Predictions of the model and implications for the gas content of globular 
clusters are discussed. 
Subject headings: globular clusters; general — stars: neutron — X-rays: general 



L INTRODUCTION 

Globular clusters contain at least two varieties of X-ray 
sources, differentiated by their luminosities. There are a dozen 
bright X-ray sources observed in the Milky Way globular clus- 
ter system, with Lx ^ 10-'''- lO-'^ erg s"' (Deutsch, Margon, & 
Anderson 2000). Nearly all of these higher luminosity sources 
show type I X-ray bursts (Lewin, van Paradijs, & Taam 1993), 
and it is, therefore, clear that they are accreting neutron stars 
(NSs) in binary systems. In fact, 7 of the 12 bright sources have 
well-measured or constrained binary periods (Deutsch, Mar- 
gon, & Anderson 2000). There is also a population of dim 
cluster X-ray sources (DCXSs), with Lx ^ 10^' - lO^^'* erg s"', 
where the lower limit is set by the detection sensitivities. To 
date, > 30 of these dim sources have been observed, primar- 
ily with the ROSAT and Einstein satellites (Hertz & Grindlay 
1983a,b; Rappaport et al. 1994; Johnston & Verbunt 1996). The 
nature of the DCXSs remains a mystery, largely due to the fact 
that sufficiently accurate positions have not been available to al- 
low for unambiguous optical or radio identifications in the host 
clusters. 

In a campaign to search for cataclysmic variables (CVs) in 
globular clusters, HST observations have discovered possible 
optical counterparts to several of the DCXSs in the clusters sur- 
veyed (Grindlay & Cool 1993; Cool et al. 1995). The absolute 
position accuracy of the HST counterparts is < 0."5, while the 
ROSAT positions of the DCXSs are no better than -- 2" (Ver- 
bunt & Hasinger 1998). In addition, the DCXSs are concen- 
trated toward the centers of their respective clusters, where the 
stellar density is high. High stellar density and source concen- 
tration, coupled with the disparity between HST and ROSAT po- 
sition accuracies, make source confusion a serious issue. Cur- 
rent and pending Chandra observations are likely to resolve the 
problem of source identification. 

As a population, the DCXSs are characterized by their low 
luminosities, temperatures kT ^ 0. 1-0.5 keV (Johnston & Ver- 
bunt 1996) for an assumed blackbody spectrum, and the fact 
that they typically lie within a few core radii from the cen- 
ters of their host clusters (Johnston, Verbunt, & Hasinger 1996; 
Verbunt & Hasinger 1998). However, this information pro- 
vides only weak constraints for models of the DCXSs. The 
hypothesis that the DCXSs are CVs (Hertz & Grindlay 1983a) 
is appealing, but the X-ray luminosities of CVs in the Galac- 



tic disk typically lie in the range Lx ^ 10 - 10 (see Fig. 8 
in Verbunt et al. 1997), falling short of the higher luminosities 
observed among the DCXSs (Verbunt, van Paradijs, & Elson 
1984). Other suggestions regarding the nature of the DCXSs 
include low-mass X-ray binaries in quiescence (Verbunt, van 
Paradijs, & Elson 1984), and rapidly rotating NSs that have not 
yet been detected as millisecond radio pulsars (see Becker & 
Triimper 1999 for X-ray luminosities of known millisecond ra- 
dio pulsars). In this work we propose an alternate model that 
may explain some fraction of the DCXSs. 

A NS moving through a gaseous medium will be able to 
accrete some of this material via the Bondi-Hoyle-Lyttleton 
(BHL) process (Hoyle & Lytdeton 1941; Bondi & Hoyle 1944; 
Bondi 1952), converting a fraction of the gravitational energy 
of the accreted gas into X-ray radiation. This physical picture 
was first proposed by Ostriker, Rees, & Silk (1970), and later 
extended by Treves & Colpi (1991) and Blaes & Madau (1993), 
for NSs in the Galactic disk, where the NS is assumed to ac- 
crete from the interstellar medium or from the gas in a giant 
molecular cloud. We propose the application of this basic idea 
to globular clusters to explain the DCXSs. This model has the 
rather minimal requirements that (i) there be a population of 
single NSs able to accrete from the intracluster gas, and (ii) the 
ambient gas density and temperature lie in the right range to 
yield the requisite accretion luminosities. 

One major point of uncertainty regarding BHL accretion by 
the Galactic disk population of single NSs is their distribution 
of space velocities (see Hansen & Phinney 1997). If the NS 
is moving supersonically through a gas, with speed v, then the 
BHL accretion luminosity is proportional to v~^, implying that 
an uncertainty in the speed distribution is considerably ampli- 
fied when computing the luminosity function of the accreting 
sources. The extension of the BHL accretion model to NSs 
in globular clusters does not suffer from this indeterminacy. 
For the more centrally concentrated globular clusters (e.g., 47 
Tuc, M15, NGC 6397) the central relaxation time is typically 
< 8 Gyr, so that at the current epoch the stars in the cluster 
core region should be nearly in thermal equilibrium (see Wat- 
ters, Joshi, & Rasio 2000) and the speed distribution for objects 
of a certain mass should be roughly Maxwellian for speeds less 
than the cluster escape speed. With a well-defined speed distri- 
bution, it is straightforward to compute the shape of the lumi- 
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nosity function for the isolated, accreting cluster NSs (lACNs). 
In this regard, a model qualitatively similar to our own (using a 
Maxwellian speed distribution), but applied to the Galactic cen- 
ter, was proposed by Zane, Turolla, & Treves (1996) to explain 
the diffuse X-ray emission from the Galactic center. 

In §2 we compute the accretion luminosity for a NS moving 
through a gaseous medium. We discuss the range of intracluster 
gas density and temperature required to explain the DCXSs. In 
§3 the luminosity function for the lACNs is derived under the 
assumptions of a Maxwellian distribution in NS speeds, and a 
constant gas density and temperature over the region of interest. 
Caveats concerning our model are discussed in §4. Finally, in 
§5 we summarize some of the fundamental points of the lACN 
model. 

2. ACCRETION FROM THE INTRACLUSTER GAS 

A NS of mass Mns, moving with relative speed v through a 
gas of ambient density p and sound speed Cs, accretes at the 
BHL rate (see Foglizzo & Ruffert 1997), 

.2^3/2 



M-4^(GM„,)Vv^ + cf)- 



(1) 



If we define V = {y^ + c^)^/^ and p = m^n, where mp is the pro- 
ton mass and n is the hydrogen number density, then for a NS of 
mass Mas = 14 Mq and radius R^s =10 km, the corresponding 
X-ray luminosity is given by 

Lx = eGM^MR-] 
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where e < 1 is the efficiency for converting gravitational energy 
into X-ray radiation, and the density is in units of cm"^. 

In order for the lACN model to successfully describe the 
DCXSs, the predicted X-ray luminosity must naturally span the 
range of the observed luminosities. This requires that the re- 
sultant speed V not be too large, and that the gas density not 
be too small. For those globular clusters where it is possible to 
measure the radial velocity dispersion, a, of the central popu- 
lation of stars, one finds cr ^ 5 -20 km s"' (Dubath, Meylan, & 
Mayor 1997). Generally, it is the light from giants, with mass 
^ 0.8 Mq, that dominates the sample used in the velocity dis- 
persion measurements. Energy equipartition implies that the 
one-dimensional velocity dispersions a^i and a„i2, for objects 
of mass m\ and m2, are related by m\a^^ = m2(J^2 (Binney & 
Tremaine 1987). If the dynamical temperature of the NS pop- 
ulation is the same as that of the turn-off stars with measured 
velocity dispersions, it follows that the one-dimensional veloc- 
ity dispersion of NSs is <Tns < 15 km s"', so that the speed it- 
self should not present a problem for the lACN model of the 
DCXSs. The density of the gas and its thermal speed are more 
uncertain, and these uncertainties are tied to the larger problem 
of the gas content in globular clusters, which we now discuss. 

A star with the cluster turn-off mass (^ 0.8 Mq) sheds most 
of its ^ 0.2 Mq envelope as it ascends the AGB, with a wind 
outflow speed of ~ 10-20 km s"' (Knapp & Morris 1985). 
Therefore, if the central escape speed of a globular cluster is 
> 20 km s"^ mass should accumulate in the central regions 
of the cluster. For a cluster that contains a mass M^ of stars 
within a few core radii from its center, the expected total rate 
of mass-loss from stars is Mwind '^ 10"^(Mc/10^ Mq) Mq yr"' 
(see Knapp, Rose, & Kerr 1973; Tayler & Wood 1975). As 
a globular cluster passes through the midplane of the Galaxy, 
it is expected that any accumulated gas will be swept out by 



ram pressure (Tayler & Wood 1975). In the r^ 10^ years be- 
tween midplane crossings, a globular cluster with M^ =10^ Mq 
could thereby accumulate of order 100 Mq of gas. This esti- 
mate is reduced considerably (by more than a factor of 10 in 
some cases) in more sophisticated treatments of the gas-flow 
problem in globular clusters which include heating and cooling 
processes (Scott & Rose 1975; Vandenberg & Faulkner 1977; 
Knapp et al. 1996). 

Direct observational searches for cluster gas in the form of 
molecular, neutral, and ionized hydrogen have yielded only 
non-detections (Smith et al. 1990; Smith, Woodsworth, & 
Hesser 1995; Knapp et al. 1996), implying upper limits to the 
total gas content in the range ^ 0.1 - 10 Mq. In a search for 
ionized hydrogen in six globular clusters, Knapp et al. (1996) 
found Mh+ ^0.1 Mq within about one core radius for the clus- 
ters observed, implying upper limits nH+ ^ 50-100 cm"-'. A 
density of 100 cm"-', however, gives a BHL accretion luminos- 
ity of lO^'* erg s"' , which is certainly high enough to explain the 
brightest of the DCXSs. 

A simplistic argument can be made to determine a lower limit 
to the density of the intracluster gas. Suppose that the wind 
from each star in the central region of a globular cluster is al- 
lowed to flow freely, unimpeded by the gravity of the cluster 
and interactions with the wind from other stars. Furthermore, 
assume that the spatial distribution of stars is uniform and con- 
fined within a sphere of radius r^. In this approximation, it can 
be shown that for r < ri. 
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where the scale r* = 0.5 pc has been chosen to represent the 
characteristic radius of a dense globular cluster core. When 
the cluster gravity and the interaction between stellar winds is 
taken into account, we suspect that the gas density can be sub- 
stantially larger than this value, although a hydrodynamic in- 
vestigation is certainly merited. 

At least two lines of argument suggest that the temperature 
of the intracluster gas should be > lO'* K, with a sound speed 
Cs - 10 (r/lO^* K)'/2 km s"^ Only one or two hot, post-AGB 
stars are required to provide sufficient ultraviolet flux to leave 
the gas near the cluster center in a warm (T ^ lO'* K), photoion- 
ized state (see Osterbrock 1989; Knapp et al. 1996). However, 
a simple calculation shows that only of order one in ten clusters 
should contain a hot, post-AGB star, owing to the short lifetime 
(< 10^ yr) of the phase. We return to this point in §4 and dis- 
cuss the consequences of a predominantly neutral intracluster 
gas for the lACN model. In addition to discrete ionizing energy 
sources within the cluster, complete thermalization of the out- 
flow speeds of ^ 10-20 km s"', via collisions between stellar 
winds, would also translate to a temperature of '^ 10"* K. 

We proceed under the assumption that there are no processes 
whereby the gas is continuously heated so that it maintains a 
very high temperature (> 10^ K), or is removed so efficiently 
as to render moot the discussion of the intracluster gas and the 
lACN model. A number of energetic gas removal mechanisms 
have been suggested (for a list, see Smith, Woodsworth, & 
Hesser 1995). Of these mechanisms, the most powerful is per- 
haps the sweeping/heating of the intracluster gas by relativis- 
tic winds or low-frequency radiation from millisecond pulsars 
(Spergel 1991). Millisecond pulsars are known to be abundant 
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in some globular clusters (see Camilo et al. 2000), but there is 
no direct evidence yet that the pulsar wind sweeping mecha- 
nism is actually operative. 

3. THE LUMINOSITY FUNCTION 

The derivation of the luminosity function (LF) is simplified 
considerably if the assumption is made that the gas density, n, 
the sound speed, c^, and the one-dimensional velocity disper- 
sion of NSs, (Tns, are constant over the region of interest. These 
assumptions should be reasonable within the core region of a 
globular cluster. Given the NS speed distribution, p(v), the X- 
ray LF is obtained by computing p{L) = \dv/dL\p(v), where we 
have dropped the subscript on Lx- We assume that the speed 
distribution in the central region of a globular cluster is well- 
represented by a Maxwellian: 



/ 2 V 
p(v)dv= J -exp{-^v^/al,)dv. 

V /I C pi; 



(4) 



This is certainly an idealization, since a more realistic speed 
distribution would vanish beyond the cluster escape speed. 
However, we would like to stress that the LF derived here is 
for illustrative purposes only. There are a number of important 
processes which may lead to a LF that is quite different from 
the one given below (see §4). 

For convenience, we define the following dimensionless 
quantities: 



£ = L/Lmax ; M = v/Cs ; H = 0-ns/Cs , 



(5) 



where L^ax = 47re(GM„s)^i?„s pc^, -^ is the maximum accretion 
luminosity (eq. [2], with v = 0). In terms of the Mach number 

M{i) = (r2/3 - 1)1/2^ the LF is 



p(e)d£ = tr 
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Fig. 1. — The luminosity function (top panel) and the cumulative distribu- 
tion (bottom panel) are shown for three different values of the characteristic 
Mach number, fi = {0.5, 1.0,2.0}. 



In Figure 1 we plot the LF and the corresponding cumulative 
distribution for three plausible values of the characteristic Mach 
number, ^. For the values of /i shown in Fig. 1, the LF peaks at 
£ > 10"^ and then drops off rather sharply for decreasing £. Fur- 
thermore, note that the LF is not well-fit by a power-law over 
more than one decade in luminosity. Therefore, if we suppose 
that the bright end of the LF is being observationally selected, 
then this model shows that a rising observed LF for decreasing 
luminosity does not necessarily mean that we are seeing the "tip 
of the iceberg," as inferred by extrapolating a power-law LF. 



4. DISCUSSION 

In deriving the LF (eq. [6]) using the standard BHL formula 
for the accretion rate (eq. [1]), we have made some rather re- 
strictive assumptions. The simple BHL formalism assumes that 
the flow is hydrodynamic (i.e., collisional), that the gas den- 
sity and temperature are uniform well beyond the accretion ra- 
dius of the NS, and that the interaction between the flow and 
the magnetic field of the NS can be neglected. Each of these 
simplifying assumptions is carried over into the derivation of 
the LF. A more realistic treatment of the accretion process is 
beyond the scope of this paper, but we can discuss certain phe- 
nomenological consequences of lifting some of the aforemen- 
tioned restrictions. 

Roughly speaking, the gas will accrete hydrodynamically 
if the mean free path of the atomic or molecular constituents 
of the gas is shorter than the BHL accretion radius (Begel- 
man 1977; Alcock & Dlarionov 1980), /?acc = 2GM„sV"^ '- 
3 X 10'"* {M/Mq)(V /IQ km s"')"^ cm, where V is the resultant 
speed in eq. (2). For neutral hydrogen, the mean free path is 
~ lO'^ n~' cm. So, if the ambient gas is predominantly neutral, 
then a density > 100 cm"^ is required for hydrodynamic accre- 
tion. On the other hand, the mean free path of a proton in an 
ionized gas is orders of magnitude shorter than the atomic mean 
free path for T ^ 10'* K (Alcock & Dlarionov 1980). Therefore, 
for globular clusters where the gas is unionized, the BHL accre- 
tion process might not be relevant. For simplicity, we restrict 
our discussion to those globular clusters where it is likely that 
the gas near the cluster center is photoionized, 47 Tuc being a 
case in point (see O'Connell et al. 1997). 

The accretion flow onto a rotating, magnetized NS will dif- 
fer substantially from the case where the NS is nonmagnetic. 
At the very least, it is expected that the flow will be inhibited in 
the former case, either due to a relativistic pulsar wind or some- 
thing akin to the "propeller" mechanism (Dlarionov & Sunyaev 
1975). It is likely that a detailed consideration of the influence 
of the magnetic field on the infalling plasma, coupled with dis- 
tributions in field strength and rotation frequency, will yield a 
broader LF than the one shown in Fig. 1. 

There are essentially two factors which determine how much 
gas is present in a globular cluster at the current epoch: (i) the 
mass and central escape speed of the cluster, and (ii) the or- 
bit of the cluster and its current position and velocity. Massive 
globular clusters with a large central escape speed will tend to 
accumulate more gas between crossings through the Galactic 
disk. Also, a cluster will likely contain more gas at the current 
epoch if it is near the top of, or on the descent from, a moder- 
ately high-altitude Galactic orbit. 

If the gas density is highly inhomogeneous, e.g., strongly 
enhanced around a discrete number of stars with large mass- 
loss rates, the LF derived in §3 would not be valid. On the 



PFAHL & RAPPAPORT 



other hand, a more clumpy gas distribution may help to ex- 
plain the presence of DCXSs at more than a few core radii 
from the center of their host cluster (see Johnston, Verbunt, 
& Hasinger 1996), where the density of an otherwise smooth 
background of gas should be markedly reduced below its cen- 
tral value (Vandenberg & Faulkner 1977). In this case, only 
when a NS passes within < 0.1 pc of a mass-losing AGB star 
might the density be sufficiently large to yield an accretion lu- 
minosity > 10""' erg s"^. 

Several of the DCXSs show evidence for significant X-ray 
variability over timescales of ^ 1 yr (Verbunt & Hasinger 1998; 
Verbunt & Johnston 2000). Within the lACN model, such time 
variability may arise due to spatial nonuniformity in the density 
or temperature of the ambient gas over length scales < 2 AU 
(for a NS speed of 10 km s"'). Hydrodynamic instabilities, pos- 
sibly induced by the NS magnetic field, may lead to variability 
on much shorter timescales. Large uncertainties associated with 
each of these mechanisms prevent us from making any clear 
predictions for the time variability of the X-ray luminosity of 
the lACNs. However, we note that numerical simulations of 
BHL accretion show a strong tendency toward nonsteady be- 
havior (see Benensohn, Lamb, & Taam 1997). In contrast, 
there is no obvious reason why an isolated rotation-powered 
pulsar should show any significant variability in the mean X- 
ray intensity (averaged over possible X-ray pulsations), which 
may detract from the millisecond pulsar hypothesis if many of 
the DCXSs are shown to be single and variable. 

5. SUMMARY AND CONCLUSIONS 

The model we have presented for the DCXSs is rather 
generic, since the requirements of the model should not be dif- 
ficult to satisfy, at least in some globular clusters. If it turns 
out that the lACN model explains a significant fraction of the 
DCXSs, then much can be learned about the population of sin- 
gle NSs in globular clusters, the properties of the intracluster 
gas, and the BHL accretion process. On the other hand, if very 
few of the DCXSs are explained by this model (perhaps be- 
cause most of the DCXSs are shown to be millisecond pulsars 
or binaries), then strong constraints can be placed on the gas 
content of globular clusters and/or the number of single NSs 
able to accrete efficiently from the intracluster gas. 

We list a number of key points regarding the BHL accretion 
scenario and the lACN population in globular clusters. 

1 . Globular clusters with large central escape speeds which 
are presently high above the Galactic plane should contain the 
largest proportion of lACNs (see §4). The globular cluster 
47 Tuc is a prime candidate in this regard. The central es- 
cape speed and height above the Galactic disk for 47 Tuc are 
^ 60 km s"' and ~ 3.2 kpc, respectively. We also note that 
there is weak observational evidence for the presence of intra- 
cluster gas in 47 Tuc, based upon diffuse UV emission from the 
central region of the cluster (O'Connell et al. 1997), as well as 
diffuse soft X-ray emission from a possible bow shock resulting 
from the interaction between the cluster gas and the low-density 
Galactic halo gas (Krockenberger & Grindlay 1995). In addi- 
tion, we note that the dispersion measures for 20 of the millisec- 
ond pulsars in 47 Tuc exhibit a range of ±0.2 cm"^ pc around 
a central value of ~ 25.5 cm"^ pc (Freire et al. 2000). Given 
that these pulsars are distributed within the inner ^ 1 pc of the 
cluster, and if we assume that the pulsars are spread along the 



line of sight by approximately this same amount, this suggests 
that the electron density within the central region of the cluster 
could be ^ 0.2 cm"-'. 

2. If the accreted gas is thermalized at the surface of the 
NS, we would expect a blackbody spectrum, with temperature 
kT - 50 (L/10^2 gjg g-iy/4y-i/4 gY^ ^^^^^ y jg jjjg fraction of 

the NS surface onto which material is accreted. For a 50 eV 
black body, the apparent B magnitude is ^^ 32 at a distance of 
2 kpc. Therefore, if the accreting NSs radiate as blackbodies, 
and if a negligible amount of the soft X-ray and UV radiation 
is reprocessed into optical light, then these objects should not 
have detectable optical counterparts. 

3. Theoretical studies indicate that the spectra of the lACNs 
should be nearly blackbody (see Alme & Wilson 1973; Zane, 
Turolla, & Treves 2000). However, based on how poorly under- 
stood are the X-ray spectra from better studied, more luminous 
accreting NSs, we would not be surprised if the X-ray spectra 
of actual BHL accreting NSs turned out to differ markedly from 
Planckian. 

4. We would like to stress that if the lACN model we pro- 
pose accounts for a very small fraction of the DCXSs, then 
rather stringent constraints may be placed on the properties 
of the intracluster gas and/or the number of NSs able to ac- 
crete. For instance, if we assume a gas temperature of 10^ K, 
then the maximum accretion luminosity of a nonmagnetic NS 
is Lx ^ IQ^^n erg s~' (eq. [2]). If no lACNs are seen above 
the detection threshold, Lx ^ 10^" erg s"', then the gas density 
must be n < I0~^ cm"^, potentially making this the most sensi- 
tive available test for the presence of intracluster gas. 

5. If accretion onto magnetized, rapidly rotating, NSs is 
strongly inhibited, then an absence of lACNs could also imply 
that recycled pulsars dominate the population of single NSs. 
This is certainly not an untenable hypothesis, since it is likely 
that a NS will undergo at least one mass and angular momen- 
tum accretion episode in a binary system over the lifetime of 
the cluster (Rasio, Pfahl, & Rappaport 2000; Pfahl, Rappaport, 
& Rasio 2000, in preparation). 

6. A strong density enhancement accompanies the large 
mass-loss rates and slow wind speeds from AGB stars. This 
suggests that some DCXSs may by spatially correlated with 
AGB stars if the lACN model is successful. Therefore, it may 
be worthwhile to calculate a projected 2D correlation function 
for DCXSs and AGB stars for those globular clusters that have 
been observed with both Chandra and HST. 

7. Finally, as an aside, we note that the M^ dependence of the 
BHL accretion luminosity (see eqs. [1] and [2]) implies that a 
black hole of mass ~ 10 Mq could be as much as lO-' times 
more luminous that an accreting NS for the same ambient gas 
density and temperature (see Grindlay 1978). There are theoret- 
ical and observational reasons to believe that there are very few 
such objects in globular clusters (see Sigurdsson & Hernquist 
1993); however, those that remain might be easily detectable 
near the cluster centers even if they are not in binary systems. 

This work was supported in part by NASA ATP grant NAG5- 
8368. We are grateful to Al Levine for a careful reading of the 
manuscript and helpful comments, and to Fred Rasio for many 
stimulating discussions about globular clusters. We also thank 
the anonymous referee for bringing point 7 in §5 to our atten- 
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